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Abstract 

We explore an extension of the Inert Doublet Model which allows also for CP violation in 
the Higgs sector. This necessitates two non-inert doublets. The lightest neutral scalar of 
the inert doublet is a candidate for dark matter. Scanning over parameters we preserve 
the abundance of the dark matter in agreement with the WMAP data. We also impose all 
relevant collider and theoretical constraints to determine the allowed parameter space for 
which both the dark matter is appropriate and CP is violated. In addition we find regions 
where the cut-off of the model originating from naturality arguments can be substantially 
lifted compared to its Standard Model value, reaching ~ 2 — 3 TeV. 
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1 Introduction 



It is widely believed that the standard model (SM) of electroweak interactions is only an 
effective low-energy theory valid below a certain energy scale A, which is supposed to be 
of the order of 1 TeV. This view is based on the fact that radiative corrections, Smi, to the 
Higgs boson mass squared {ml), tend to increase the mass up to A, implying a necessary 
fine tuning. This is the celebrated little hierarchy problem or the "LEP paradox" [T]. In 
order to retain a meaningful perturbative expansion above ~ 1 TeV, a high level of fine 
tuning between mh and A is necessary to suppress 6m\ relative to m\, see e.g. [2]. Other 
well know problems of the SM are the lack of a candidate for the dark matter (DM) and 
too little CP violation (CPV) that could make the electroweak baryogenesis viable [3]. 
In that context the SM scalar sector has to be modifed as the phase transition within a 
single Higgs doublet is too slow for baryogenesis [3]. 

Our goal here is to outline a model that ameliorates the little hierarchy problem 
(lifting A at least to ~ 2 TeV) while providing extra sources of CP violation needed for 
baryogenesis as well as a realistic abundance of dark matter. We will focus on extending 
the Higgs sector of the SM by adding extra Higgs doublets since that could also help to 
make the electroweak phase transition fast enough, see 

In general there are two possibilities to alleviate the little hierarchy problem: (i) 
suppression of the radiative corrections, and/or (ii) an increase of the Higgs mass. Well 
known examples of the first strategy are supersymmetric extensions of the SM, however 
in fact that could be achieved to some extent even through very modest means, e.g. by 
introducing just one extra real scalar singlet to the SM [3] (although some fine tuning 
is necessary). The second possibility has recently been followed by Barbieri, Hall and 
Rychkov [6]H The idea is to introduce an extra scalar doublet i] (the inert doublet) 
which does not couple to the SM fermions as a consequence of a Z2 symmetry: 1] —>■ —r] 
(all other fields are neutral under the symmetry). Since it is required that (ry) = 0, the 
symmetry remains unbroken, therefore it can provide a DM candidate. Since the inert 
doublet contributes to gauge-boson two-point Green's functions the SM-like Higgs boson 
could be as heavy as 400 — 600 GeV, ameliorating the little hierarchy problem this way. 
Also the DM constraints could be satisfied choosing masses of the scalar and pseudoscalar 
components of the inert Higgs of the order of 80 GeV. The analysis of |9] reveals also 
another solution for the DM candidate, such that scalar and pseudoscalar masses are much 
heavier, > 500 GeV. The model in very simple terms avoids the little hierarchy problem 
for the lighter (inert) scalars as they just do not couple to fermions, so in particular not to 
the top quark. The only drawback is that the model is so restricted by the Z2 symmetry 
that it does not allow for CP violation in the Higgs potential, this is the issue that we 
would like to address here. There are two simple extensions of the Inert Doublet Model 
(IDM) that can accommodate CP violation: 

• Combining the standard 2HDM with an inert scalar doublet 77 (in other words 
replacing the SM Higgs doublet of the IDM by two doublets). 

• Adding a complex singlet scalar to the IDM. 

This paper is devoted to the first of these extensions. It is worth emphasizing that this 
scenario is not just a simple sum of the 2HDM and the IDM. Although some theoretical 

^The model was proposed earlier in [7] as a possible solution to the DM problem, its collider phe- 
nomenology was then discussed in [Sj . 
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and experimental constraints which are apphcable here, are (to leading order) identical to 
those of the 2HDM (just because the inert doublet does not couple to fermions), there are 
also important exceptions. These concern the oblique parameters T and S, the amount 
of DM and the issue of positivity of the scalar potential (vacuum stability): 

• The extra inert degrees of freedom transform as an SU(2) doublet, so they couple 
to the vector bosons and therefore contribute to the oblique parameters, modifying 
the standard 2HDM predictions. 

• The neutral components of the inert doublet are candidates for DM, and since the 
inert doublet couples to the 2HDM doublets, the amplitudes for DM annihilation are 
in general influenced in a non-trivial way by the extension of the non-inert sector. 

• Even under simplifying assumptions, the scalar potential for the 2HDM extended 
by the inert doublet has a rich structure, so that the condition for positivity is much 
more involved than a simple superposition of conditions needed for the 2HDM and 
the IDM separately. 

The paper is organized as follows. In Sec. [2] we introduce the model and define some 
notation. In Sec. [3] we define some benchmarks for the inert sector, and in Sec. H] we 
present the strategy adopted to search for allowed regions in the parameter space of the 
model. Sees. E] and are devoted to reviews of theoretical and experimental constraints. 
In Sec. [7] we show some regions of parameters of the model that are compatible with all 
the constraints, and in Sec. Owe summarize. 

Technical details on positivity, CP conservation and necessary basis transformations 
are collected in appendices A, B and C. 



2 Inert-plus-two-doublet model: IDM2 
2.1 The potential 

Introduction of two doublets, $1^2 leads in general to Flavor-Changing Neutral Currents 
in Yukawa couplings. To avoid those one can impose an extra Z2 symmetry such that 
$1 — ^ — $1 and ur —ur (all other fields are neutral). The model then has Z2 x Zg, 
where the first factor is the inert-doublet Z2: f] ^ —rj (all other fields are neutral). The 
potential reads 



y($i, $2, r]) = Vi2($i, $2) + V3{t]) + Vusi'^i, $2, v) 



(2.1) 



where 



^2(^1, $2) = ~ + ml2<l>l<^2 + [m?2$l$2 + h.c] } 



+ ^($I$l)^ + y ($^$2)' + A3($I$l)($5$2) + A4($I$2)($5$l) 



Ai 

r 

1 

+ 2 



A5($I$2)' + h.c. 



A 



(2.2) 
(2.3) 
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Fl23($l,$2,^?) = Aii33($l$l)(r/^??) + A2233(*2*2)(r?M 

+ A 

1331 

{^[r]){r]^^^) + A2332($k)(^^$2) 



+1 



Al313('^'k)' + h-C. +i A2323($k)'+h.C. 



(2.4) 



Here, A1133, A2233, A1331 and A2332 are real, whereas A1313 and A2323 can be complex. 
Disregarding $2? the correspondence with the notation of [U] would be 



(m^, A^, A1133, A1331, A1313) ^ (/i2, 2A2, A3, A4, A5). 



(2.5) 



In Vi($i, $2) we have allowed for soft breaking of Z2 in order to preserve the chance 
of CP violation in the potential while we do not allow for any breaking of Z2 in order to 
have a stable lightest component of 77 as a DM candidate. Note that, as a consequence 
of the unbroken Z2, there is no mixing in mass terms between $12 and 77. It is worth 
realizing that, since 77 does not couple to quarks, there are no constraints on the charged 
inert Higgs mass from the 6 ^ 57 decay. 



2.2 2HDM mass eigenstates 

In the (non-inert) 2HDM sector of the model, we denote the doublets (in a basis where 
both have a v.e.v.) 

( (^l + ^l + %)/v^)' "^'^ ( (^^2 + ^2+%)/v^)' ^^-^^ 

and adopt the mixing matrix R, defined by 



HA 






H2 




V2 






Vl3 



satisfying 



RM'^R' = A<Lg = diag(M^, M|, M|), (2.^ 



(2.7) 



and parametrized in terms of three rotation angles as 

Ci C2 Si C2 S2 

R= \ -{C1S2S3 + S1C3) C1C3-S1S2S3 C2S3 I (2.9) 

-Ci S2 C3 + Si S3 — (Ci S3 + Si S2 C3) C2 C3^ 

with Ci = cosaj, Sj = sinaj. In Eq. f l2.7p . 773 = —smj3xi + cos/5x2 is the combination of 
Xi's which is orthogonal to the neutral Nambu-Goldstone boson. Here, tan/5 = V2/V1. 

We also define a mass parameter /i^ = {v'^ /2viV2)^eml2, and note the following useful 
relation: 

Imm^2 = ImA5fif2. (2-10) 
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2.3 Inert-sector mass eigenstates 

Components of the inert doublet are defined as follows 



(5 + L)/v^)- (^-^^^ 
The masses of the inert scalars will be given by expressions analogous to those of [6l [9] : 

Ml = + i [Al + Al- Al,,,] , (2.12) 
where we have introduced the abbreviations 

^ch = All33f? + A2233 ^^2) 
Aq = Ai33it;i + A2332 ^^2? 

A%,^ = Re Ai3i3 vf + Re A2323 (2.13) 
Adopting the simplifying assumptions (denoted "dark democracy" ) 

Aa = A1133 = A22335 
Afe = A1331 = A2332 5 

Ac = Ai3i3 = A2323 (real), (2.14) 



the masses can be written as: 



Ml = ml + i(Aa + \b + = + |(Ab + AJ^;^ 

Ml = + \{Xa + = M,'± + i(Afe - X^y. (2.15) 

As a consequence of the assumptions (12.141) . there are no trilinear couplings H^tj^S or 



2.4 Stability of the potential 

The condition for positivity of V is discussed in Appendix A for the general potential, 
Eq. ( 12.11) . In our numerical applications we will limit ourselves to the case of "dark 
democracy" defined in (12.141) . We find that for this special case, Eqs. ( 1A.16P and ( 1A.32I) 
must be satisfied for positivity. However, we restrict ourselves even further by requiring 
V12, V3 and V123 separately to be positive. Then, in addition to the familiar constraint on 
^2 in [ini [H] and V3, we obtain the following condition: 

A, > max(0,-2Ab,-Ab±Ae), (2.16) 

implying < M^±. This amounts to a strong constraint on the splitting of the inert- 
sector spectrum, not present in the full treatment of positivity. 

The input parameters in the inert sector are defined by specifying scalar masses 
(M5,Ma,Mj^±) together with m^, so that the quartic couplings Aa, A{, and Ac can be 
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determined via (12.151) . Here we will consider cases (the profile 3 is the only exception, see 
Sec. 13.21) . with masses ordered as follows: Ms < M4 < M^±. This is motivated by the 
fact that a positive contribution to the electroweak precision observable T from the inert 
sector (see Sec. E]) makes it easier to allow for heavy 2HDM masses. For this case it is easy 
to show that the ordering and the positivity condition (I2.16P leave a certain non-empty 
allowed region in the space of (Aa,Afe, Ac), namely Aa > 0, A^ < 0, \a > 2|Afe| together 
with Ac < and |Ab| > |Ac|. Then the requirement of having the right amount of dark 
matter (see next section) imposes additional constraints on the masses (or equivalently 
on the quartic couplings), resulting in a relatively small region of allowed (Aa, A;,, Ac). If 
Ms < M^± < Ma, then the constraints are less tight. 

3 Benchmarks 

In numerical studies we will assume that 5* is the lightest neutral scalar, Ms < M4. The 
first profiles have the DM candidate around 75 GeV, a favoured value P, |9]. For the 
heavier neutral partner, we consider a few options at, or slightly above Ma = 110 GeV, 
which is the lower limit compatible with LEP2 data |12j . 

In a recent study by Lopez Honorez et al. [S], the splitting amongst inert-sector scalar 
masses was kept fixed, while a scan over the parameter fi2 (corresponding to our m^) and 
the DM particle mass {Ms in our notation) was performed. In the region of heavy DM 
particles, a detailed study has also been performed in [13j. 

We shall instead consider two discrete sets of "dark" profiles, to be specified below. 
In the first set of "dark profiles", we keep the DM particle light. For trial values of the 
masses and (inspired by the results of [U]), we estimate the amount of dark matter 
from micrOMEGAs2.2 [T^ITS] . in the IDM version developed by Lopez Honorez et al. [9]. 
That version has only one Z2-even doublet, whereas we here consider two such doublets, 
$1 and $2, (with many more "free" parameters). Therefore the calculation of amplitudes 
for various DM annihilation channels is more complicated. 

The mass of the charged partner, M^± , and the inert-sector mass parameter have 
been chosen such that a reasonable amount of dark matter is obtained for at least one set 
of the parameters (13. 6p . In the context of dark matter, the parameter is important, 
since a particular choice for the inert neutral scalar masses Ms and Ma together with 
constrain (Aq, A;,, Ac) which, in turn, are responsible for the annihilation of dark matter 
into the visible sector, see [9]. In the original IDM, the essential parameter determining 
the trilinear coupling among two DM-particles and the SM Higgs boson, is 

Az. = |(Aa + A5 + Ac). (3.1) 

In the present model, the trilinear coupling between two S'-particles and a neutral 
Higgs boson is determined by 

AL(fi?7i + V2r]2)SS. (3.2) 
Projecting out the coupling to a particular neutral Higgs boson, we find: 

SSHj : Fssj^L, with Fssj = cos PRji + sin (3 Rj2, (3.3) 

where the pre-factor satisfies \Fssj\ < 1, since R is unitary. In particular, Fssi = cos(/3 — 
ai) cos 02- 
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Similarly, the four-point coupling involving two S'-particles and two neutral Higgs 
bosons is determined by 

^[(A, + A, + KM + r,l) + (A, + A, - Ae)(x? + xl)]SS. (3.4) 

Projecting out the couphng to the hghtest Higgs boson, we find: 

SSH^H^ : i(Ai - KRl,) = - Aesin^as). (3.5) 

In our estimates of the DM density, we shall follow the approach of [9j, taking A^ as the 
relevant parameter. The DM values presented in this work are based on the full tree- 
level calculation performed using micrOMEGAs2.2 with all intermediate and final states 
originating from the rich 2HDM structure of the modelJl 

3.1 2HDM masses 

For each of the dark profiles, we consider the following mass parameters of the (non-inert) 
2HDM sector: 

Set A : (Ml, Ms) = (100, 300) GeV, /i = 200 GeV, (3.6a) 
Set B : (Ml, Ms) = (200, 400) GeV, ^^ = 400 GeV, (3.6b) 
Set C : (Ml, Ma) = (400, 500) GeV, /i = 400 GeV. (3.6c) 

A non-zero value for fi is adopted, in order to accommodate the unitarity constraints lim- 
iting quartic couplings. We avoid degeneracy of Mi and M2, since that would be a source 
of potential difficulties to produce CP violation, see Sec. 15.31 for a detailed discussion. 

3.2 Light DM particle 

We shall consider the following "light" DM profiles: 



Profile 


1 : 


Ms 


= 75 GeV, 


Ma = 


110 


GeV, 


M^± 


= 112 GeV, 


Profile 


1' : 


Ms 


= 77 GeV, 


Ma = 


110 


GeV, 


M^± 


= 112 GeV, 


Profile 


2 : 


Ms 


= 75 GeV, 


Ma = 


120 


GeV, 


M^± 


= 125 GeV, 


Profile 


3 : 


Ms 


= 75 GeV, 


Ma = 


120 


GeV, 


M^± 


= 85 GeV, 


Profile 


4 : 


Ms 


= 100 GeV, 


Ma = 


110 


GeV, 


M^± 


= 115 GeV, 


Profile 


5 : 


Ms 


= 120 GeV, 


Ma = 


125 


GeV, 


M,± 


= 130 GeV. 



(3.7) 



3.3 Heavier DM particle 

We also consider some profiles where the two neutral inert-particle masses are higher, and 
rather close, another domain favoured by [U]: 

Profile 11 : M^ = 500 GeV, M^ = 501 GeV, M^± = 502 GeV, 

^Within the approximate treatment of positivity, and keeping only the lowest mass state of the 2HDM, 
the full result can be mimicked by an appropriate tuning of m^, which in turn amounts to a tuning of 

Aa. 
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Profile 12 : Ms = 600 GeV, Ma = 601 GeV, M^± = 602 GeV, 
Profile 13 : Ms = 800 GeV, = 802 GeV, M^± = 804 GeV, 

Profile 14 : = 1000 GeV, Ma = 1002 GeV, M^± = 1005 GeV. (3.8) 



The latter profiles have a high degree of degeneracy among the masses. This is required 
in order to have the correct amount of dark matter, but will also minimize the contribution 
to the electroweak observable T. For a detailed study of the IDM in the high-mass regime, 
see [13]. 

4 Search Strategy 

In general, our goal will be to verify that within the IDM2 model one can accommodate 
the following features 

• a "heavy" lightest neutral scalar (so that the naturalness problem is alleviated), 

• at least one neutral scalar odd under Z2 consistent with the present limits on the 
DM abundance, 

• CP violation in the potential (introduced via -^s)- 

First we choose the following input parameters for the 2HDM: tan/5. Mi, M2, Mh± 
and II, together with (ai, 0:2, 03). All remaining parameters of the 2HDM sector are 
calculable in terms of those chosen above [TU] (see also [I^). For the inert sector we 
choose M5 ~ 70 - 80 GeV or > 500 GeV (DM candidate). Ma, M^± and (needed for 
determination of relevant quartic couplings in V123 which are necessary for calculating the 
DM abundance). 

Then the following strategy will be applied while determining allowed regions in the 
parameter space of the model: 

• We fix tan j3, Mh±, sets of 2HDM masses (Mi, M2, /i) and profiles of the inert-sector 
parameters ("dark profiles", Ms, Ma, M^±). 

• Next we check if for a given choice of dark profile (M5, Ma, M^±) and (Mi, M2, jj), 
there exists a value of such that the model predicts the right order of magnitude 
for the dark matter abundance. 

• Then we scan over the mixing angles (ai, 02, cts) of the 2HDM, requiring that: 

— The naturalness is alleviated for each Higgs boson 



where M denotes a generic Higgs boson mass, 5M^ stands for the top-quark 
contribution to the one- loop correction to M^, at is a calculable coefficient in 
terms of the mixing angles etc. The cut-off should be chosen to be in the TeV 
region, e.g. modestly A ~ 2 TeV. The fine-tuning parameter D is to be chosen 
according to our aesthetic standards. 



\5M^ 
M2 




(4.1) 
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— Based on the experience gained from [6] we restrict the scan to heavy (> 
100 GeV) Higgs bosons in the 2HDM sector (so that the httle hierarchy problem 
could be that way reduced, that would be an analog of the heavy SM Higgs 
of [6]). It is worth noting that here some possible tension between parameters 
emerges. It may appear as we increase Higgs masses in the 2HDM sector trying 
to retain small quartic constants. Even though the masses could be raised by 
increasing /i^ ~ the potential, nevertheless the mass of one scalar would 
still remain ~ v (see e.g. [IE])- In order to increase its mass some combination 
of quartic couplings in Vu will have to be large, therefore checking the unitarity 
in the 2HDM sector is essential to guarantee that Aj's remain in a perturbative 
regime. 

— Remaining experimental constraints are satisfied. 

— CP is violated, i.e. aj are far enough from their CPC limits. We use the 
invariants Jj, i = 1,2,3 [19] as a measure of CP violation. Eventually we plot 
an average and/or maximum (with respect to the mixing angles a,) for |Im Ji| 
and the electron electric dipole moment (EDM) to illustrate the strength of 
CP violation. 



5 Theoretical constraints 
5.1 The little hierarchy 

As an order-of-magnitude estimate for radiative corrections to neutral Higgs boson masses, 
we consider the contributions that arise from top-quark loops: 

^^l = -4^^'(«? + S') for ^■ = 1'2,3 (5.1) 
where aj and aj are defined in (3.21) of [20] : 

ttj = aj = — 5.2) 

Similarly, for the charged Higgs particles we find 

5Ml^ = -^,^l^ooi'(3. (5.3) 

Since the inert doublet does not couple to fermions there is no hierarchy problem for S 
and A. 

We adopt the following simple condition 



J ll < D \ U±l < (5 4) 

M] ' 

with the amount of fine tuning parametrized by D. For the resulting cut-off we will chose 

A = min(Aj, Ah±). (5.5) 

This quantity A will in general be most constrained by the value of Mi (unless when 
a\ -|- a\ is small). For an alleviation of the hierarchy problem, we would like to have A 
"large" compared to the Higgs masses. 
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5.2 Perturbativity and Unitarity 

To preserve perturbativity we shall impose the following conditions on quartic and Yukawa 
couplings of neutral and charged Higgs bosons 

\, \aj\, a,- , cot Aa, Af,, Ac < 47r. (5.6) 

Since we will consider tan/3 > 0.5, the last three conditions will always be satisfied. We 
also impose unitarity on the Higgs-Higgs scattering amplitudes 



5.3 CP Violation 

Since our intention here is to outline a model which would possess CP violation in the 
Higgs potential we shall discuss this issue in more detail. The magnitude of CP violation 
can be quantified in terms of the invariants introduced by Lavoura and Silva [THj. However 
here we prefer to adopt the more general, basis-independent approach of Gunion and 
Haber and calculate the invariants Ji, J2 and J3 of [M]. They state (Theorem 4) that 
the Higgs sector is CP-conserving if and only if all J, are real. The calculations of these 
quantities are straightforward, we end up with the following result, valid for our choice of 
basis: 

2 2 

ImJi = -!^(Ai-A2)ImA5, (5.7) 

ImJ, = [{(X, - X, - - \X,\^) vt + 2{\^ - \,)Re X.vfvl 

- ((A2 - A3 - A4)' - I Asp) vf\ Im A5, (5.8) 

2 2 

ImJa = ^(Ai-A2)(Ai + A2 + 2A4)ImA5. (5.9) 

We note that since we have chosen a basis with real v.e.v.'s, there is no CP violation 
when ImAs = (it should be realized that Imm^g ImAs are not independent here, 
see (12.101) ). Then an interesting question arises: Is it possible, for ImAs 7^ to have no 
CP violation? It turns out that the answer is "yes" , as will be discussed in the following. 

The simultaneous vanishing of the three Im Jj implies that CP is conserved. This can 
happen for five distinct cases: 



• Case A: ImAs 7^ with vi = (tan/3 00). 

• Case B: Im A5 7^ with f 2 = (tan /3 = 0). 

• Case C: ImAs = 0. This corresponds to Ai^ (see (12. 8p ) being block diagonal, and 
the diagonalization performed in terms of only one rotation angle. 

• Case D: ImAs 7^ with Ai = A2 and vi = V2- 

• Case E: ImAs 7^ with Ai = A2, vi 7^ V2 and (Ai — A3 — A4)^ = |Asp. 
Some comments are here in order (details are discussed in Appendix B): 
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• When all three masses are degenerate, CP is conserved because ImAs = 0. This 
corresponds to Case C. 

• When there is only partial mass degeneracy, Mi = M2 < M3 or Mi < M2 = M3, 
there are instances of CP conservation corresponding to cases C, D and E. 

• There are also instances of CP conservation in the mass non-degenerate case Mi < 
M2 < M3 corresponding to cases C, D and E. 

The above discussion shows that, in terms of our input parameters (tan/3. Mi, M2, Mh± 
and /i, together with (01,02,03)) there exist various non-trivial locations such that CP 
is conserved even though ImAs 7^ (the case of ImAs = is relatively obvious). Since 
our intention is to build a model that allows for a substantial amount of CP violation, 
we would like to show regions of parameter space where that indeed happens. However, 
in light of the above discussion, the determination of such locations can not easily be 
performed analytically. 

The chance for successful electroweak baryogenesis is the crucial motivation for our 
discussion of CP violation. However, without a dedicated analysis of baryogenesis (which 
is beyond the scope of this project) it is hard to estimate the amount of CP violation that 
is necessary. Therefore, we have adopted the following strategy to illustrate the strength 
of CP violation which is available in the model. We plot both the electron electric dipole 
moment and the invariant |Im Ji| (which are physical quantities) in the region of the 
parameter space allowed by all the other constraints. In order to estimate the amount of 
potential CP violation we show both averaged and maximal values of and |Im Ji| (the 
choice of Ji (as opposed to J2 and J3) is arbitrary, however we recall that it is sufficient for 
CP violation to have just one of the Ji complex). Large splitting between averaged and 
maximal values indicates the potential for CP violation hidden in the appropriate choice of 
mixing angles Oj. Of course, in a realistic situation (having the prediction for electroweak 
baryogenesis within the model) we would need to have d^ below the experimental upper 
limit and nevertheless enough CP violation for successful baryogenesis. 

In Appendix B, in Figs. [7] and El we show how allowed regions in the (oi, 02, as) space 
are distributed. 

6 Experimental constraints 

We here review various experimental constraints that will be imposed on the model. 

T and S: We adopt the results from [25], [26j in order to calculate T and S within our 
model. Since for the Higgs fields we will use a basis in which only $1 has non-zero vev (so 
called Higgs basis) some necessary transformations must be performed, see Appendix C 
for details. For the model discussed here the rotation matrix O defined by Eq. (59) of 
I25] reads 




O 



I On O12 Oi3 \ 

O21 O22 O23 

O31 O32 O33 

1 

\ 1 / 



(6.1) 
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and the mixing matrix V of that paper becomes 



V 



On 
O21 + iOu 




O12 
O22 + iOz2 




Oi3 
O23 + ^033 






(6.2) 



with U = 1. Therefore, U'^V = V, and 



( 1 



\ 



-lOn 
1 

-Z0i3 
%0x2 






^Oi3 
1 

-lOn 





-20i3 

-iOi2 

1 










1 







i 

1 / 



(6.3) 



Since this is block diagonal, the contribution to T (and to S") from the inert doublet is 
additive (the inert fields must always appear in pairs, there is no interference between the 
"visible" and the inert sector at the one- loop order): 



T — ^2HDM + 



—-4- ^ [F(M2±, M|) + F(M2±, M\) - F(Ml, Ml 

167r sm B-wrriy^ ' ' 



0] (6.4) 



where Q;e.m.r2HDM = Ap2HDM IS given by (63) of [25]. In our case, the matrices ( 1C.3P and 
( ]C.6p should be adopted. Similarly, S can be obtained from the results given in |26j. We 
impose the bounds |Ar| < 0.10, |AS| < 0.10 [27J, at the \-a level. 

We note that T, which is our main concern, gets a positive contribution from a split- 
ting between the masses of charged and neutral Higgs bosons, whereas a pair of neutral 
ones gives a negative contribution. In fact, since the function F is symmetric in its two 
arguments, these two opposite-sign contributions cancel in the limit when the charged 
boson is degenerate with either of the two neutral ones. 

-Bo — -Bo mixing: Due to the possibility of charged-Higgs exchange, in addition to 
exchange, the -Bq — -Bq mixing constraint excludes low values of tan (5 and low values of 
Mh± [281 EHl EOj. Here we follow the procedure of [31]. 

-B — > Xs'y'- The b ^ s'j transition may also proceed via charged-Higgs exchange, 
so some regions of the parameter space with low values of tan /? and Mh± are excluded. 
The exact region of exclusion is sensitive to higher-order QCD effects [321 ESI EH ES] , and 
roughly excludes Mjj± < 300 GeV. Again, we follow the approach of ^T\. 

B — > tUt-X: The charged Higgs contribution may substantially modify the branching 
ratio for B — >■ rVrX [36] . The measurement [37i| of B{B —>■ rz/^-X) leads to the following 
constraint 



tan/5 



< 0.53 GeV" 



(6.5) 



at 95% CL. This is in fact a very weak constraint. A more recent measurement gives 
B{B^ — s> tVt-) = (1.79±0.71) X 10~^ [3H], where we have added in quadrature symmetrized 
statistical and systematic errors. With the SM prediction of (1.59 ± 0.40) x 10""^, 



exp 



B{B- 
B{B- 



1.13 ±0.53. 



(6.6) 
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Within the framework of the 2HDM, one finds [391 

ml 



^^H 2HDM 



Then the data imply that two sectors at large values of tan (3 and low values of M^± are 
excluded. 

B — > DtUt- Measurements HQ] of the ratio 



pexp B{B Dru^) 

can also be used to constrain the coupling of the charged Higgs to the r, more precisely 
ta.nP/MH±. It thus restricts large values of tan/3 and low values of Mh± [H], in a region 
of parameter space similar to the one following from B rUrX, but is considerably 
stronger. 

LEP2 non-discovery: The non-discovery of a Higgs boson at LEP2 imposes a bound 
on how strongly the lightest one can couple to the Z and to bb. Useful results for this 
constraint are available in table 27 of |i42j. Adopting the standard notation 

(^z{H^bb) = ^zfn^bb) ^ '^z{H^bb) (6-9) 
one can approximately parametrize the upper limit on Cz{H^bb) fro^ the table as follows: 

r 0.05 for 12 GeV < < 80 GeV, 
Cl(^H^f,l) < < 0.1 for 80 GeV < < 90 GeV, (6.10) 
[ 0.2 for 90 GeV < Ma < 110 GeV. 

For e+e" — > ZHi the coefficient C'^z{H^bb) given by (4.3) of [20j: 

C\h.^u>) = {c^Rn + s^Ruf \ {Rli + slRl,) . (6.11) 

Then the constraint fl6.10p limits the parameter space. 

Ri,: The branching ratio for Z ^ bb is also affected by Higgs exchange. As noticed in 
|20] . the contributions from neutral Higgs bosons to R^ are negligible, however, charged 
Higgs boson contributions, as given by ^3], Eq. (4.2), exclude low values of tan/? and low 
Mh±. Experimentally Rb = Tz^bb/'^z^h^d = 0.21629 ±0.00066 [27]. It is easy to see that 
the correction STz^bh implies the following change for Rb. 

6Rb=^^^{l-Rb) (6.12) 

J- Z^had 



where Tz^had = (1.7444±0.0020) GeV [27J. Since 6Tz^bb is known within the 2HDM 
so is 6Rb. We require 

5Rb < 0.00066, (6.13) 

corresponding to STz^bb = 1-47 MeV at the 1-a level. 

Muon anomalous magnetic moment: Since here we are considering heavy Higgs 
bosons {Mi > 100 GeV) therefore, according to [SI [31] , the 2HDM contribution to the 
muon anomalous magnetic moment is negligible even for tan /? as large as ~ 40. 
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Electron electric dipole moment: The bounds on electric dipole moments con- 
strain the allowed amount of CP violation of the model. We adopt the bound P3| (see 
also [46J): 

141 < 1 X 10~2^[ecm], (6.14) 

at the 1-0" level. The contribution due to neutral Higgs exchangelf], via the two-loop Barr- 
Zee effect [IE], is given by Eq. (3.2) of [16] in terms of the neutral-sector mixing matrix 
O, defined in [IH], and related to our R via 

Oil Oi2 Oi-A I -R33 -R23 -RisX 

O21 O22 023 = Rzi R21 Rii . (6.15) 
O31 O32 033/ \Rm R22 -R12/ 

This rotation matrix O should not be confused with the one appearing in Eq. fl6.ip . 

7 Results 

Subject to the limitations discussed above, we may now scan over the parameter space, 
imposing the constraint 

X' = l]x?< 5.99, 95%C.L., (7.1) 

as appropriate for identifying allowed regions in two dimensions, (tan/?, Mj^±). The sum 
runs over all the experimental constraints discussed above. 

Our results will be given in terms of contour plots of various quantities of interest. 
Regions corresponding to values being confined within certain intervals are indicated by 
a colour coding as indicated. The external contour shows the maximal region consistent 
with the experimental and theoretical constraints we adopt (i.e., without the naturality 
condition imposed on A/-\/D, unless explicitly stated). 

The model discussed here contains many parameters (masses, mixing angles, etc.). 
In projecting down our results to a lower- dimensional space, we have decided to favour 
the more "physical" parameters tan (3 and Mh± together with neutral scalar masses (and 
/i). In this section, most of the plots show (for fixed inert masses, M5,M^,M^±, fixed 
2HDM neutral masses, Mi,M2 and /i) allowed regions in the (tan/5, M/^±) space. The 
remaining parameters, the neutral- Higgs-sector mixing angles, ai, ^2, and 03, which are 
not specified in those plots, have been averaged over, or a maximum has been extracted. 
Thus, for each allowed point in (tan/3,M//±) there exist a's such that all constraints are 
satisfied. 

7.1 Light Dark Matter particle 

We first consider the set of "dark profiles" for which the DM particle is light, at ~ 75 GeV, 
see Eq. f l3.7l) . These profiles may be compared with the parameters considered by Lopez 
Honorez et al. [9], in their Fig. 5, where it was shown that the correct amount of dark 
matter may be obtained with dark-matter mass of the order of 50-80 GeV, for splittings 
of the order 10-50 GeV. 

•^Neglecting the CKM mixing, the charged Higgs contribution to d,, vanishes in the 2HDM (with softly 
broken Z2) up to two loops ^47j. 
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Inert masses : 75, 1 1 0, 1 1 2 GeV ^2'^- ^ 00. 300, 200 GeV 




Inert masses: 75, 110, 112 GeV ^2'^'- ""OO, 300, 200 GeV 




10 tanp 1 10 tanp 



Figure 1: Top panels: Relative amount of CP violation. Bottom panels: Electron electric 
dipole moment in units [e 10~^^ cm]. Left: average; Right: maximum value. Inert-sector 
masses: (75, 110, 112) GeV; 2HDM-sector masses: (Mi, M2, /i) = (100, 300, 200) GeV. 

In Table [H we summarize the results for the DM density Qh"^ for the different combi- 
nations of dark profiles of Eq. (13. 7p . and non-inert-sector parameters of Eq. (13.61) . These 
are obtained from running micrOMEGAs f!M [TH] . as described in Sec. El For the quartic 
couplings in the inert sector (denoted here by A^/2), we use the value 0.1 (the amount 
of dark matter is practically independent of this coupling [9]). Here, "OK" means that 
a value within 3a of the WMAP value 0.1131 ± 0.0034 [50], can be found for a suitable 
choice of m^. For these profiles of light DM particles, typically values ~ 30 — 50 GeV 
are required. 

Profile 1' corresponds to the allowed horizontal band in Fig. 9 of Lundstrom et al. [12], 
in the sense that we find solutions for Sets B and C, corresponding to heavier non-inert 
Higgs particles. This band is rather narrow: for Ms = 75 GeV we find too much dark 
matter, at Ms = 79 GeV too little. It appears shifted by a couple of GeV, with respect to 
the results of [12], presumably due to the use of a different DM code, micrOMEGAs [HI [15] 
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DPI 


DPI' 


DPS 


DP4, DP5 




(75,110,112) 


(77,110,112) 


(75,110,85) 


(100,110,115) 


A (100,300,200) 


OK 


< 0.07 


< 0.09 


< 0.01 


B (200,400,400) 


> 0.24 


OK 


OK 


< 0.01 


C (400,500,400) 


> 0.24 


0.08-0.09 


< 0.10 


< 0.01 



Table 1: Dark-matter density nh^ for different "dark profiles", DPI, DPI', DPS, DP4, 
DP5 with mass parameters {Ms, M^, M^±), and different mass-parameter sets A-C with 
mass parameters (Mi, M2,/i), all in GeV, see Sec. [Sl 

vs. DarkSusy [5T]. 

All these profiles (except Profile S) have the charged particle heavier than both the 
neutral ones, in order that the inert sector makes a positive contribution to T (see Sec. E]). 
This, in turn, allows for a higher value of the "ordinary" neutral Higgs particle [H] that 
enters the electroweak fits. 

Profiles 4 and 5 give Qh"^ < 0.01 (consistent with Figs. 8 and 9 of [12]) and will not 
be considered any further. 

Dark profile 1. We start discussing the case (M5,M^,M^±) = (75,110,112) GeV. 
Among the three sets of mass parameters considered for the visible sector, only Set A 
(Ml, Ms, /i) = (100, SOO, 200) GeV gives a reasonable value of nh"^. With of the order 
45-55 GeV we obtain Qh"^ ~ 0.10 — 0. IS. Higher values of yield too high values of Qh'^. 
For Profile 1 and Sets B and C, the value is too high, nh'^ > 0.240 The allowed region 
is shown in Fig. [H The figure is obtained by scanning over tan/3 from 0.5 to 50, and 
over Mh± from SOO GeV (barely above the B Xg'-f cut-off) to 700 GeV. For each point 
in the {tan /3, Mh±) plane, a scan over mixing angles {ai, 02,0^3) is performed, analyzing 
all models compliant with the constraints, showing contours of the CP-violating quantity 
|Im Ji| and the electron electric dipole moment (EDM). These two quantities provide two 
different measures of the amount of CP violation. The outer contours delineate the regions 
within which consistent solutions are found for one or more sets of mixing angles Oj. The 
averages presented in the left panels are obtained by averaging over all sets of (ai, 02, c^s) 
for which the experimental constraints are satisfied according to Eq. (17. ip . Similarly, the 
maxima in the right panels correspond to maxima of absolute values, obtained from the 
same scans over allowed sets of (ai, ^2, as). 

In the upper part of the figure we display |ImJi|, while the lower one shows the 
electron electric dipole moment, \de\- In the left panels averages over allowed sets of 
are displayed, while the right ones show extrema. There is no obvious correlation between 
I del and |Im Ji|, other than both having maxima in the interior of the allowed region. 

In Appendix B we discuss the correlations of angles for which viable solutions are 
found. These are shown separately for "small tan/3" and "large tan/3". In general, it 
is easier to accommodate CP violation at low values of tan /3 than at higher values. For 
high values of tan the allowed values of these parameters tend to accumulate near the 
limits where H2 is odd under CP (0^2 ~ 0, as ~ ±7r/2). 

For the parameters considered here, (Mi, M2, /i) = (100, SOO, 200) GeV, the viable 
models are constrained to tan/3 ~ 0.5 — 6. The cut-off at "high" tan/3 is mostly due to 

''For Sets B and C the lightest 2HDM Higgs is relatively heavy so that SS annihilation via an inter- 
mediate Hi is suppressed, and too much dark matter would survive. 
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the unitarity constraint on the neutral-Higgs sector. We also note a cut-off at Mh± ~ 
650 GeV. This, on the other hand, is due to the electroweak constraints, T in particular. 

The fact that the average electric dipole moment is rather small compared to the 
maximum value, means that large parts of the tan (3-Mh± space would remain viable 
even if the experimental constraint on \de\ should become significantly tightened. 

Inert: 75, 110, 112GeV M,. f^2'^' 100,300, 200 GeV 

4000 
3500 
3000 
2500 
2000 
1500 
1000 
500 

iJWU 

1 10 tanp 

Figure 2: Maximum reach in A/\^ [GeV]. Inert-sector masses: (75,110,112) GeV; 
2HDM-sector masses: (Mi, Ms, /i) = (100, 300, 200) GeV. 

The little hierarchy. In this Fig. [1], no constraint from the fine-tuning consideration 
in Sec. 15.11 is imposed. In Fig. [2] we plot the quantity A/a/D of Eq. (15.51) . which by 
definition is minimized over the three neutral Higgs particles and the charged one (for 
fixed mixing angles ««). This is next maximized over the ai (for fixed tan/3 and Mh±) 
before being plotted in this figure. For the cases considered, we note that A/v^ reaches 
up to around 3 TeV. The dominant reason why this varies with tan/3 and Mii± is that 
the factor 

= a] + ~a] (7.2) 

of Eq. (15. ip . for which solutions are allowed, varies. Specifically, when the factor is 
small, k/\/D will be larger. This happens when sin ai cos 0^2 and cos/? sin 0^2 are both 
small (cf. Appendix B). 

Accordingly, imposing the constraint (15. 4p . with K/\/D = 2 TeV, we obtain the more 
restricted allowed regions shown in Fig. [SI Compared with Fig. [T] (where there is no 
constraint on A), we see a dramatic reduction of the allowed parameter space. For Mi = 
100 GeV, only a small region of Mjy± ~ 300 — 350 GeV and tan /? ~ 2 — 6 survives. 
Naturally, this is dominantly caused by the condition (I5.4p applied to Mi. 

Dark profile 1'. This profile has a slightly higher value of M^, and thus gives allowed 
models also out to higher values of the non-inert Higgs mass. Mi, corresponding to the 
"horizontal band" in Fig. 9 of Ref . . The results for Sets B and C are shown in Fig. HI 
As compared with profile 1, the heavier S here allows more annihilation to (off-shell) WW 
pairs, and thus an acceptable value of Vth? can be found. 

For Set B, there are solutions over most of the explored tan (3-Mh± plane, for Mh± ~ 
400 - 500 GeV and M2 = /i = 400 GeV they reach all the way out to tan (3 ~ 50. 
For Set C, on the other hand, the situation is different. Apart from a region around 
tan/3 ~ 1 — 2, solutions are only found within two widely separated bands in the tan/3- 
Mh± plane. There is one band at Mh± ^ 400 GeV and another at Mh± < 550 — 650 GeV, 
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Inert masses: 75, 110, 112 GeV ^2'^- lOO. 300' 200 GeV 




10 tanp 1 10 tanp 



Figure 3: Allowed regions for A/\^ = 2 TeV (outer contours) and electron electric dipole 
moment in units [e 10^^^ cm]. Left: average; Right: maximum value. Inert-sector masses: 
(75, 110, 112) GeV; 2HDM-sector masses: (Mi, M2, /i) = (100, 300, 200) GeV. 



with a region of no solutions in-between. This is clearly due to the interplay or partial 
cancellation of contributions to T, positive from the charged-Higgs boson differing in mass 
from some neutral ones, and negative from pairs of neutral ones having different masses. 

Dark profile 2. In this case the allowed models are located in essentially the same 
part of the tan (3-Mh± plane as for "dark profile 1", with similar values for |Im Ji| and 
the electron electric dipole moment. 

Dark profile 3. This profile differs from the previous ones in having the charged 
boson of the inert sector lighter than the heavier neutral one. With this modification 
from Profile 1, we find acceptable solutions also for heavier Higgs bosons in the non-inert 
sector (Sets B and C). The allowed regions in the tan (3-Mh± plane are quite similar to 
those shown in Fig. Hlfor Profile 1', and also the values of the electric dipole moment are 
similar. 

The little hierarchy. If we here impose the constraint 05.41) . with A/^/D = 2 TeV, 
there is a considerable reduction in the allowed parameter space for Set B, whereas for 
Set C essentially the whole parameter space survives, as shown in Fig. [5] for Profile 1'. 

7.2 Heavier Dark Matter particle 

The second set of "dark profiles" has a heavier DM particle. For this case, it was shown 
(see Fig. 6 of [9]) that the amount of dark matter can satisfy the WMAP constraint 
[50] provided the splitting in the inert sector is small. We confirm that finding and list 
in Table [2] approximate values of flh'^ obtained from micrOMEGAs [HI [T3] with suitable 
choices of m^. The values vary quite a bit with the choice of this parameter. As mentioned 
above, the small mass splitting in this inert sector also leads to a small (hence acceptable) 
modification of T. 

Dark profiles 11, 12, 13 and 14. For all these cases, we find solutions for all the 
three considered sets of non- inert mass parameters of Eq. (13. 6p . The allowed regions in the 



18 



Inert masses: 77, 110, 112 GeV 



M ^: 200, 400, 400 GeV 




10 tanp 

Inert masses: 77, 110, 112 GeV 



Ml, M ^i: 400, 500, 400 GeV 




tanp 



tan|3 



Figure 4: Electron electric dipole moment in units [e 10 cm]. Left: average; Right: 
maximum value. Inert-sector masses: (77, 110, 112) GeV; 2HDM-sector masses: (Mi, M2, 
fi) = (200, 400, 400) GeV and (400, 500, 400) GeV. 



tan /3-Mh± plane and the corresponding values for the electron electric dipole moment are 
very similar to those displayed in Figs. [Hand HI For Set B, with (Mi, M2, /i) = (200, 400, 
400) GeV, the allowed region extends to large values of Mh± without coming into conflict 
with the EW precision data, and to large values of tan/3 without coming into conflict 
with unitarity. However, the large-tan j3 part is constrained to Mh± ~ 400 — 500 GeV by 
the electroweak precision data, T in particular. This region roughly corresponds to the 
decoupling limit. For Set C, there are two disconnected, allowed regions. 

The little hierarchy. We show in Fig. El the resulting ranges in A/-\/zJ, as defined by 
Eq. (15.51) . This is seen to reach well beyond 3 TeV, in particular for Set B. Comparing sets 
of the 2HDM mass parameters it is clear that, in agreement with our expectations, the 
maximal value of the cut-off A/ \/D grows with the mass scale of the 2HDM model. This 
illustrates our strategy to ameliorate the little hierarachy problem by lifting the lowest 
visible Higgs mass. As already mentioned, the allowed regions in this case are very similar 
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Inert masses: 77, 110, 112 GeV ^2'^- 200, 400, 400 GeV 




Inert masses: 77, 110, 112 GeV ^2'^- 400, 500, 400 GeV 




10 tanp 1 1° tanp 



Figure 5: Allowed regions for A/^/D = 2 TeV (outer contours) and electron electric dipole 
moment in units [e 10~^^ cm]. Left: average; Right: maximum value. Inert-sector masses: 
(77, 110, 112) GeV; 2HDM-sector masses: (Mi, M2, /i) = (200, 400, 400) GeV and (400, 
500, 400) GeV. 

to those corresponding to light dark profiles shown in Figs. [1] and HI This is consistent 
with the fact that the inert sector influences the experimental constraints only through 
T which is sensitive to inert scalar-mass splitting. Since the splitting is of the same order 
(although it is larger for the light dark profiles), therefore, it is not surprising that the 
allowed regions are similar. Consequently, if we impose a cut on K/\/D = 2 TeV we find 
allowed regions very similar to those found for the lighter dark profiles, see Figs. [3] and [5l 
In general, the heavier non-inert Higgs states affect the amount of dark matter. How- 
ever, within the range of parameters explored, their effect can be compensated by a 
retuning of the soft mass parameter m^, which for fixed masses yields a calibration of the 
trilinear inert-non-inert Higgs coupling A^. 
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DPll 


DP12 


DP13 


DP14 




(500,501,502) 


(600,601,602) 


(800,802,804) 


(1000,1002,1005) 


A (100,300,200) 


~ 0.09 


^ 0.10 


^ 0.09 


~ 0.10 


B (200,400,400) 


~ 0.09 


^ 0.10 


^ 0.09 


~ 0.10 


C (400,500,400) 


~ 0.09 


~ 0.10 


~ 0.09 


~ 0.10 



Table 2: DM density nh^ for different "dark profiles", DPll, DP12, DP13, DP14 
with mass parameters {Ms, Ma, M^±) and different 2HDM mass-parameter sets A-C, 
(Ml, M2, /i), all in GeV, see Sec. O 

8 Summary 

We have explored an extension of the Inert Doublet Model [SI El IE] (IDM), made by 
replacing the SM Higgs doublet sector by a 2-Higgs-Doublet Model (2HDM) in order to 
accommodate CP violation in interactions of neutral Higgs bosons. This model has four 
inert-sector scalars: two neutral, and a pair of charged ones, as well as three "ordinary" 
neutral scalars and an accompanying pair of charged ones. The latter five are those of 
the familiar 2HDM. Our motivation here was not only to have CP violation in the scalar 
potential, but also to provide a candidate for dark matter and to ameliorate the little 
hierarchy problem by lifting Higgs boson masses (thereby increasing the cut-off A). 

We have estimated the amount of dark matter that is predicted by the model, adopt- 
ing the code micrOMEGAs [HI [IS]) checking all relevant theoretical and experimental con- 
straints. Solutions were found both for "light" DM particles, with a mass around 75 GeV, 
as well as for heavier ones, of the order of a few hundred GeV. In both cases, the splitting 
between the masses of the charged and the neutral inert scalars must be small in order 
to reproduce the right amount of dark matter. In the case of heavier dark matter the 
splitting is tiny implying nearly vanishing contribution to the T parameter. 

As we have noted at the very end of Sec. 17.21 the regions allowed by the condition 
A/^/D > 2 TeV for lighter and heavier dark matter profiles are similar. This is a conse- 
quence of a very small contribution to T from the inert sector (because of the small mass 
splitting). Since the inert sector infiuences the experimental constraints only through 
T, the allowed regions are similar. This observation illustrates an important difference 
between this model and the original IDM, where the charged inert scalar has to be con- 
siderably heavier than the neutral inert scalars in order to provide a contribution to T 
that can compensate a large negative SM contribution (for a heavy SM Higgs boson). 
Here the role of the inert sector is mainly restricted to providing a candidate for the 
dark matter, while T can be made consistent with the data utilizing only the freedom of 
the 2HDM sector with negligible contribution from the inert sector. That freedom was 
not available in the original inert model. In our approach the little hierarchy problem is 
softened by increasing the Higgs boson masses, this is possible mainly through the 2HDM 
sector alone. 

In general, and in agreement with our expectations, the naturality arguments favour 
heavier 2HDM masses. For instance, if one requires > 2 TeV, then as seen in the 

top panel of Fig. [HI only a small region with tan /5 ~ 2 — 5 and Mfj± ~ 320 — 380 GeV 
is allowed for (Mi, M2, fx) = (100,300,200) GeV. In contrast, for heavier 2HDM masses 
(the middle panel in the same figure) a much larger region remains with tan /? ~ 1 — 40 
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Inert: 500, 501, 502GeV M,. f^^'^' 100,300, 200 GeV 
700r 




tanp 

Inert: 500, 501 , 502 GeV M„ M^, n: 200,400,400 GeV 
700r 



> 

V 

a 




Inert: 500, 501 , 502 GeV M„ M^, \i: 400,500,400 GeV 

700r 




Figure 6: Maximum reach in A/vD [GeV]. Inert-sector masses: (500,501,502) GeV; 
2HDM-sector masses: (Mi, M2, fi) = (100, 300, 200) GeV, (200, 400, 400) GeV and (400, 
500, 400) GeV. 



and Mh± ~ 400 - 640 GeV. On the other hand, for (Mi,M2,/i) = (400,500,400) GeV 
(the lower panel) we observe again a large allowed region, which however in this case is 
split into two disconnected regions corresponding to light {Mh± ~ 300 — 450 GeV) and 
heavy {Mh± ~ 550 — 650 GeV) charged Higgs bosons. The splitting is related to partial 
cancellation between contributions to T from charged and neutral 2HDM scalars. It is 
worth noting, that the pictures which emerge here are similar to the result obtained within 
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the ordinary 2HDM [3T] (see also [52]). As we have aheady mentioned, that similarity 
follows from the small splitting in the inert doublet-sector masses. 

For the light-DM scenario discussed in Sec. 17.1^ the experimental implications (direct 
detection and production at the LHC) were addressed in P IS, |9]. For the case of heavier 
DM particles, the prospects are more dim, since the production cross section will be rather 
small. For a detailed discussion, see [53l [Ml [I3j . 

As for the CP invariance, we observe that in general, CP can be substantially violated 
in regions allowed by the experimental and theoretical constraints. For instance, it is seen 
that the electron electric dipole moment can easily exceed the allowed value of 10"^'' [e cm] 
for appropriate choices of the mixing angles. We also note that it is easier to accommodate 
CP violation at low values of tan/3 than at higher values. For high values of tan/3, the 
allowed values of (01,02,03) tend to accumulate near the limits where H2 is odd under 
CP (02 ^ 0, 03 ~ ±7r/2). 

It is also worth realizing that the model we discuss here bears some similarity to the 
Weinberg model [55] of CP violation with natural absence of flavour-changing Yukawa 
couplings. Both models invoke three doublets, one of which has no Yukawa couplings 
while the two others couple to fermions such that no flavour- changing Yukawa couplings 
emerge. However, the important difference is that here (in order to guarantee stability of 
the dark matter candidate) we impose an extra Z2 symmetry (77 — —rj) which remains 
unbroken after spontaneous symmetry breaking since (r/) = 0. As a consequence, there is 
no mixing in the mass matrix between the charged components of $1,2 and r], therefore 
eventually r/^ is a mass eigenstate with no Yukawa couplings at all. Then CP violation 
in the charged scalar sector is the same as in a pure 2HDM with Yukawa couplings 
parameterized by the CKM matrix alone. 

It should be noted that restricting our study to the range of inert-model parameters 
satisfying (12.161) . rather than utihzing the full range defined by Eqs. ( ]A.16|) and ( ]A.32[) . 
we are clearly not able to find all allowed domains in the parameter space. However, a 
full investigation is technically much more involved and therefore computationally more 
challenging. In order to determine all the quartic couplings in the 2DHM sector (V12) 
one has to specify {Mi, M2, fi), tan/?, Mh± and also the angles (01,02,03)- To fix the 
quartic couplings between the 2HDM and the inert sector contained in V123 one has to 
know also the inert mass parameters [Ms, Ma, M^±) and m^. Then for each point in the 
parameter space (including angles over which we scan) both Eqs. flA.16p and flA.32p must 
be checked. In the presence of the large number of free parameters, that makes the analysis 
much more complicated and time consuming. Also, within the general strategy outlined 
here, the calculation of the DM abundance would be much more complicated, since for 
that purpose all the parameters must be simultaneously known. That is indeed necessary 
as cross sections for dark matter annihilation depend on the mixing angles (01,02,03), 
masses of the 2HDM scalars Mi, M2, M3 and, of course, also on the inert mass parameters 
{Ms, Ma, Mjj±) and m^. A more complete investigation adopting the general positivity 
conditions Eqs. ( 1A.16I) and ( 1A.32I) together with a more precise calculation of the dark 
matter abundance will be attempted, and reported on elsewhere. 
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A Positivity of the potential 

In order to ensure vacuum stability, the potential should be positive for large values of 
the fields, $i, $2 and rj. We study the general potential 

m,$2,r/) = ^($I$i)^ + ^($^$2)^ + ^(r/V 

+ A3($I$l)($^$2) + A4($I<f 2)($5$l) + I [A5($I$2)' + h.C 
+ Aii33($I$l)(r/"^r?) + X2233{H^2){ri^ri) 

+A 

1331 

{^\r]){r]^<!>i) + A2332($k)(^^$2) 



+ 7 



+ 



A2323($k)^ + h.C. 



|m^i$J<l>i +m^2'^5<l>2 + m^2'^|<l>2 + h.C. | 



We start by rewriting the Higgs doublets as: 

$1 = ||<l>l||<l>l, $2 = ||$2||<l>2, 



(A.l) 



(A.2) 



where and ||?7|| are the norms of the spinors, and and r) are unit spinors. By 

SU (2) invariance, only the following combinations of fields may appear: 



2^2 



i$2ir, v^v 



$5$! = ||$i|| ■ ||<l'2|| ■ l-l) , $I$2 = H^lY, 

7^t(|,^ = W^iW ■ ||?7|| ^57"^ ■ <lij , ^{l] = [t]''^i]*, 

V^^2 = \ \^2\ \ - Ml ■ ^2^ , = [V^^2]*- 

We let the norms of Eq. (\A.2\} be parametrized as follows: 

||$i|| = r cos7sin^, ||$2|| = ''"sin7sin^, Hr^H = r cos^^. 



(A.3) 



(A.4) 



The complex product between two different unit spinors will be a complex number with 
modulus less than or equal to unity, i.e. 



$5 ■ $1 = Pie 



fl' ■ <l>2 



i4>3 



(A.5) 
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Using this parametrization, we can write: 

= cos^ 7 sin^ $2*^2 = ''"^ sin^ 7sin^ rj^ rj = r"^ cos^ 6 , 
$2$! = r^cos7sin7sin^^pie*'^\ $|$2 = COS7 sin7 sin^ 6'pie"*'^\ 
rj^^i = cos 7 sin 9 cos 6p2e^'^'^ , ^\r] = r'^ cos 7 sin cos 6p2e~^'^'^ , 
7yt<|)2 = sin 7 sin 6* cos Op^e^'^^ , ^Irj = sin 7 sin ^ cos 6pse~^'^^ (A. 6) 

where r > 0, 7 e [0,7r/2], ^ e [0,7r/2], e [0, 1] and (pi e [0,27r). 
The potential can now be written as 

^ = rV4 + rV2, (A.7) 

with only the quartic, V4, part relevant for positivity: 

Vi = XiAi + A2A2 + X^As + A3A4 + A4A5 

+^1133^6 + A2233^7 + A1331A8 + A2332^9 

+Re A5 Aio + Im X5A11 

+Re A1313A12 + Im A1313A13 + Re A2323^i4 + Im A2323^i5, (A. 8) 

where 

Ai = ^ cos** 7 sin^ ^, ^2 = ^ sin'' 7 sin^ ^, 743 = ^cos''^, = cos^ 7 sin^ 7 sin^ ^, 

^5 = Pi cos^ 7 sin^ 7 sin^ 6*, = cos^ 7 sin^ 6* cos^ = sin^ 7 sin^ 6* cos^ 6', 

As = p\ cos^ 7 sin^ 6 cos^ 6^ Aq = p\ sin^ 7 sin^ 6 cos^ 6^ 
^10 = Pi cos(20i) cos^ 7 sin^ 7 sin'' 6, An = pi sin(20i) cos^ 7 sin^ 7 sin^ 9, 
A12 = P2 cos(202) cos^ 7 sin^ 9 cos^ ^, A13 = P2 sin(2(/)2) cos^ 7 sin^ 9 cos^ ^, 
Ai4 = pI cos(203) sin^ 7 sin^ 9 cos^ 6', A15 = pi sin(203) sin'^ 7 sin'^ 9 cos^ 6'. (A. 9) 

The quartic part of the potential can now be written: 

4 • 4 /) I -^2 . 4 . 4 ^ A^ 4 

= — COS 7 sm + — sm 7 sm 9 + — cos d' 

Zi Zi Zj 

+A3 cos^ 7 sin^ 7 sin'^ + A1133 cos^ 7 sin^ ^ cos^ 9 + A2233 sin^ 7 sin^ 9 cos^ ^ 

+Pi [A4 + Re A5 cos(20i) + Im A5 sin(20i)] cos^ 7 sin^ 7 sin'' 9 

+P2 [-^1331 + Re A1313 cos(202) + Im A1313 sin(202)] cos^ 7 sin^ 9 cos^ 9 

+pI [A2332 + Re A2323 cos(203) + Im A2323 sin(203)] sin^ 7 sin^ 6* cos^ 9. (A. 10) 

We minimize this expression with respect to 0i to arrive at 

f> '^l 4 • 4 /I I . 4 . 4 /I , '^Tj 4/1 

V4 — cos 7 sm p + sm 7 sm 9 -\- cos ^ 

Zi Zi Zi 

+A3 cos^ 7 sin^ 7 sin^ ^ + A1133 cos^ 7 sin^ ^ cos^ 9 + A2233 sin^ 7 sin^ ^ cos^ ^ 
+Pi ( A4 — I A5 1 ) cos^ 7 sin^ 7 sin^ ^ 
+P2 (A1331 - IA1313I) cos^ 7sin^ ^ cos^ ^ 

+P3 (A2332 - IA2323I) sin^ 7sin^ 6lcos^ 9. (A.ll) 
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Further, we minimize this expression with respect to pi to arrive at: 

Va — COS 7 sm 6 ^ sm 7 sm U ^ cos U 

2 ' 2 ' 2 



where 



+ cos^ 7 sin^ 7 sin^ 6* + Ay cos^ 7 sin^ 6^ cos^ 6* + A^ sin^ 7 sin^ 6* cos^ 6. (A. 12) 



A, = A3 + min(0,A4- IA5I) (A.13) 
\ = A1133 + min(0,Ai33i - IA1313I) (A-M) 
Az = A2233 + min(0, A2332 - IA2323I) (A.15) 

For the positivity condition to be satisfied, V4 must be positive for all combinations of 
7 G [0,7r/2] and 6 G [0,7r/2]. This is both a necessary and a sufficient condition. 



A.l Boundary points 

Some points from the parameter space give us some rather simple positivity conditions. 
We now turn our attention towards these special points. 

^ = or ^ = 7r/2 or 7 = or 7 = 7r/2 

First we consider the boundary points in the (7, 6) plane. 

The last three of these expressions have the same form as an expression already studied 
in the 2HDM [20]. Using a result from there, we end up with the following conditions: 

Ai > 0, A2 > 0, A^ > 0, X^> -VA1A2, Xy > -v^AiA^, A^ > -^/X^. (A.16) 



^\ -^1 4 , -^2 . 4 , , 2-2 

= -) = ycos 7 + ysm 7 + Aa;COS 7sm 7, 

= 0) = y sin^^ + ^cos^e + Aj^sin^^cos^e, 

= -) = ^sin^^ + ^cos^^ + A.sin^^cos^^. 



A. 2 Interior points 

What remains is to demand that V4 > also in the interior of the (7, 6) plane. Thus, 

Xy cos^ 7 + A2 sin^ 7 ^ ~o I — 2^ ~'~ ^"^1 '"^^^ 7 + A2 sin'' 7 + 2Aa; cos^ 7 sin^ 7) tan^ 1 

2 I tan u J 

Maximizing the right hand side of this inequality with respect to 6, we find that the 
maximum occurs at tan^ 6 = a/ A^/ (Ai cos^ 7 + A2 sin'' 7 + 2Aa; cos^ 7 sin^ 7) . Substituting 
this back we arrive at: 

Xy cos^ 7 + Az sin^ 1 ^ ~\J cos^ 7 + A2 sin'^ 7 + 2Aa; cos^ 7 sin^ 7) (A. 17) 
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We need to solve this inequality subject to the constraints given in (]A.16I) . Let us distin- 
guish between four different cases: 

Case a) > and > 



The inequality ( 1A.17I) is trivially satisfied. 



Case b) Aj, > and A^ < 

The left-hand side of (IA.17I) can be both positive and negative. See details in section lAlSl 
Case c) Xy < and A^ > 

The left-hand side of f lA.17p can be both positive and negative. See details in section 1X31 
Should be similar to case b) with Xy and A2 interchanged. 

Case d) (A^, < A A^ < 0) or (Aj, < A A^ < 0) 

We can now square both sides of the inequahty and reverse the inequahty sign to get 
Xy cos^ 7 + 2XyXz cos^ 7 sin^ 1 + sin^ 7 < A,,(Ai cos"^ 7 + A2 sin^ 7 + 2Xx cos^ 7 sin^ 7) 
and finally 

(A^Ai - A,^) cos^ 7 + (A^A2 - A^) sin^ 7 + 2(A^Aa; - A^A^) cos^ 7 sin^ 7 > 
which is positive definite if 



A^A, - A,A, > -y/(A,Ai-A2)(A^A2-A2) (A.18) 



A. 3 A detailed study of positivity in cases b) and c) 

We consider the inequality flA.17p : 



Aj^cos 7 + A2sin 7 > — y A^(Ai cos^ 7 + A2 sin 7 + 2Aj: cos^ 7 sin 7). (A. 19) 
Let us introduce x = tan^ 7 and transform the inequality into 



Aj, + X,x > -^JX^{Xl + A2X2 + 2A^x). (A.20) 

which must be satisfied for x > 0. In order to analyze this we will study the solutions of 
the equation 



Ay + X^x = -^X^{Xi + A2X2 + 2A^x). (A.21) 
Possible solutions of this equation are given by 



-b ± _ 4^^ 
xi,2 = ^ (A.22) 
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where a = Xr,X2 — ^l, b = 2(Ar,Aa; — Aj^A^) and c = A,,Ai — Xy. These possible solutions are 
obtained by squaring flA.2ip . However, in doing this we may introduce false solutions, 
that is solutions of 



Xy + XzX = +\ A^(Ai + X2x'^ + 2A^.x). 



(A.23) 



In order for ( ]A.20[) to be satisfied we must demand that ( ]A.21[) does not have any positive 
real- valued solutions when subject to the constraints already obtained in ( ]A.16|) . That is: 



62 - 4ac > A [(xi > A Ay + X^Xi < 0) V (xs > A Aj^ + A^Xa < 0)] 



(A.24) 



where the overline (bar) denotes negation. This expression must be true in order for 
positivity to be satisfied. The possible solutions read: 

Case b) A^ > and A^ < 

In this case a > 0, while c can be both positive, zero or negative. We have to distin- 
guish between these cases. The results are: 



cA^ 

OOA (a>-^V6> -2y/acj 



Vc < 



(A.25) 



equivalent to 

A^Ai >Xl^ [X2XI > XiXl V A^A, - XyX, > -Y^(A^Ai- A2)(A,,A2- A2) 
Case c) Aj, < and A^ > 



V A^Ai < X'y 



(A.26) 



In this case c > 0, while a can be both positive, zero or negative. We need to distinguish 
between these cases. The results are: 



aA^ 

a> A [ c> b> -2y/ac 



Va < 



(A.27) 



equivalent to 



A,A2 > A^ A AiA^ > X2X' V A,A, - A,A, > - J(A,Ai - A2)(A,A2 - A^) 



V A^As < Xl 



(A.28) 



Combining the results from cases a) - d) subject to the constraints of flA.16p there are 
remarkable(!) simplifications. We end up with 



AiA, + V^Ay>0 V A^A, - Aj,A, > (A^Ai - A2)(A^,A2 - A^) 



(A.29) 



as an additional constraint to the ones listed in (lA.16p . or expressed in a more symmetric 
form: 



X^X^ + ^/XlX, + ^/X2Xy>0 V A^A^ + AiA^ + A2A2-A^AiA2-2A^.Aj,A, < (A.30) 
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A. 4 The dark democracy 

In the so-called dark democracy (12 ■141) . = and flA.17p simplifies to 



Aj^ > — Y A^(Ai cos"' 7 + A2 sin 7 + 2Aa; cos^ 7 sin 7) (A. 31) 

which is satisfied whenever 

A, > V (a,A, - AJ > -y/(A,Ai-A2)(A,A2-A2)) (A.32) 



This constraint combined with the inequalities from ( ]A.16I) forms a necessary and sufficient 
condition to guarantee positivity of the potential. 

B CP conservation 

In the general case, the mixing of the different weak neutral states is described by the 
angles {ai, a,-, as} present in the rotation matrix R. There are three simple limits of no 
CP violation, namely when either i^i, or is odd under CP. These limits can all be 
defined in terms of a, and 0:3, as illustrated in Fig. 1 of [T7j: H\ is odd when a, = ±7r/2, 
H2 is odd when a, = and a, = ±vr/2, whereas is odd when a, = as = 0. Away 
from these limits, the model in general violates CP. (Exceptions will be discussed below.) 

We display in Figs. [7] and M how the allowed regions of the a parameter space are 
distributed for different choices of the mass parameters and different cuts on tan/3. These 
plots are obtained by collecting all points for which solutions (in terms of tan f3, Mh± , ai, a,, 03) 
are found during a particular scan (indicated by the parameters at the top), and binning 
them in the variables. In these figures, we also indicate (yellow plane and green lines) 
limits where there is no CP violation. These figures show that in general, it is much easier 
to accommodate CP violation at low values of tan f3 than at higher values. 

In addition to the CP-conserving limits mentioned above there exist also other regions 
in the parameter space which imply CP invariance. In the following we will identify those 
regions in terms of the scalar masses, tan/5 and ctj. In this process we will make repeated 
use of the following identities that follow from the orthogonality of R: 



Ml, 


= rUmI-mI 


) + rUmI-mI) 


+ Ml 


Ml, 


= Rl,{Ml-Ml 


) + RUMl-Ml) 


+ Ml 


Ml, 


= Rl,{Ml-Ml 


) + rUmI-mI) 


+ Ml 


Ml, 


= RnRuiM^- 


Ml) + R,,R„{Ml 


-Ml) 


Ml, 


= RnRniMl- 


Ml) + R,,R„{Ml 


-Ml) 


Ml, 


= Ri,Ri,{Ml- 


Ml) + R„R,,{Ml 


-Ml) 



In addition we order the masses so that Mi < M, < M,. 
B.l ImA5 = 

Whenever Im A5 = we know that CP is conserved. To see when this happens, we note 
that 

- ^s^Im As = Ml, = RiiRn{Ml - Ml) + R,iR,,{Ml - Ml) 
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Inert masses: 75, 110, 112 GeV M^, Mj, ^: 100, 300, 200 GeV 




Figure 7: Populated regions of a-space. Left: tan (3 < 1; Right: tan [3 > 3. Inert-sector 
masses: (75,110,112) GeV; 2HDM-sector masses: (Mi, Ms, n) = (100, 300, 200) GeV. 
Green lines at ^2 = (and as = 7r/2 or 0) show limits of no CP violation, with either H2 
or H3 being odd under CP. The limit 02 = vr/2 (yellow plane) indicates where Hi is odd. 
There is a corresponding one at a2 = — 7r/2 (not shown). 

Inert masses: 77, 110, 112 GeV M^, Mj, ^: 400, 500, 400 GeV 




Figure 8: Populated regions of a-space. Similar to Fig. [7] for inert-sector masses: 
(77, 110, 112) GeV; 2HDM-sector masses: (Mi, M2, /i) = (400, 500, 400) GeV. 

-^CfslmX, = Ml, = RuRisiM^ - Ml) + RsMM^ - Ml). (B.l) 

Both TVI 13 and AI23 niust vanish in order to get Im A5 = 0. (If only one of these quantities 
vanishes while the other is nonzero, that would correspond to/3 = 0or/3 = 7r/2.) That 
happens in the following cases 

• Full mass degeneracy Mi = M2 = M3. 

• Lower mass degeneracy Mi = M2 < M3, when one of the following conditions is 
satisfied 
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— a2 = vr/2ora3 = 7r/2 (This corresponds to R33 = 0.) 

— a2 = and as = 0. (This corresponds to -R31 = and -R32 = 0.) 

• Higher mass degeneracy Mi < M2 = M3, when one of the following conditions is 
satisfied 

— a2 = (This corresponds to -R13 = 0.) 

— ^2 = 7r/2 (This corresponds to Ru = and R12 = 0.) 

• Nondegenerate masses Mi < M2 < M3, when one of the following conditions is 
satisfied 

— a2 = vr/2 (This corresponds to -Ru = and R12 = and -R33 = 0.) 

— 0^2 = and a-s = (This corresponds to Ris = and -R31 = and R32 = 0.) 

— a2 = and 0:3 = n/2 (This corresponds to R13 = and i?33 = 0.) 

B.2 Partially degenerate masses 

We have seen in the previous section that full mass degeneracy implies Im A5 = 0. How- 
ever, when ImAs 7^ 0, partial mass degeneracy is still possible in some cases. We will 
point out those cases here. We begin by noting that 

MI3 - tsuaPMls = 

or, equivalent ly 

RisiRn - i?i2tan/5)(M2 - M|) + R^^iRsi - i?32 tan /3) (M^ - M|) = 0. (B.2) 

• Lower mass degeneracy Mi = M2 < M3 is allowed when i?3i — i?32 tan (3 = 0. 
(Excluding i?3i = i?32 = which would imply Im A5 = as explained in the previous 
sub-section.) 

• Higher mass degeneracy Mi < M2 = M3 is allowed when Ru — R12 tan [3 = 0. 
(Excluding Ru = R12 = which would imply Im A5 = as explained in the previous 
sub-section.) 

B.3 Ai = A2 and tan/3 = 1 

We know that even when Im A5 7^ we can have cases of CP conservation. One such case 
is when Ai = A2 and tan/? = 1. 

[Ai-A2U,=i = ^,{M\,-Ml2) 

= ^^[{Rl,-Rl^{Ml-Ml) + {Rl,-Rl^){Ml-Ml)\ (B.3) 



This expression must be zero subject to the constraint (1B.2P in order to have CP conser 



vation. Thus, when ImAs 7^ and tan/3 = 1 we have CP conservation in the following 
cases 
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• Lower mass degeneracy Mi = M2 < M3 when R^i — R32 7^ 0. 

• Higher mass degeneracy Mi < M2 = M3 when i?ii = R12 7^ 0. 

• Nondegenerate masses Mi < M2 < M3 in one of the following cases: 

— When _Rii = R12 and i?3i = i?32. (Excluding the cases which would imply 
ImA5 = 0.) This corresponds to ai —tt/A, — and 0:2 arbitrary (but not 
or 7r/2). 

— When i?ii = i?i2 and -R33 = (this implies -R31 = —R32)- (Excluding the cases 
which would imply ImAs = 0.) This corresponds to ai = 7r/4, = tt/2 and 
a2 arbitrary (but not or 7r/2). 

— When i?3i = i?32 and i?i3 = (this implies i?ii = — -Ri2)- (Excluding the cases 
which would imply ImAs = 0.) This corresponds to ai — — 7r/4, a2 — and 
as arbitrary (but not or 7r/2). 



B.4 Ai = A2 and (Ai - A3 - X^f = 1^5? 

Finally, when ImAs 7^ we also have CP conservation when Ai = A2 and (Ai — A3 — A4)^ = 
lAsp. In order to see what this corresponds to, we start by solving the equation Ai = A2 
for the parameter u. We find 



u — 



- CIMI2 



when Ai = A2. Furthermore, we find that 



(B.4) 



(Ai - A3 - A4)' 



lAsr = 









Ml 



2 



12 



-Mis 



By substituting the expression for v which is valid when Ai = A2, we arrive at 
[(Ai-A3-A4f-|A5n,^^,^ 

= ^4^2(^2 _i) + 4)(-M22 - M\^{Ml2 + - 2Mls) 

+2i^(l + t%){Mls - M\2)M\2 - 2H{\ + t%){M\s - Ml,)Ml2 

+(ti - i)(i + tinMi2r - 4^(4 - 1) [(Mi,r + iMi,r] } 

After substituting the mass matrix elements into this expression we get: 



[(Ai-A3-A4r-|A5n,^^,^ 
_ (M| - M|)2(i?3i - tpR,2) 



{(1 + 4)(i?31 + tpR32){Rz2 - 
+2tpRl^ [(i?3l(t? - Mf,) + i?32(l - 3i|)] } 
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(Mf - MinRn - t,R^ ^ ^ ^^^^^^^^^^ _ ^^^^^y 



vHia-tv 



+2tpRl^ [{Rii{tl - 3tp) + i?i2(l - 3tl)] } 

2{Mf - M^){Ml - Ml'' 



I) 

-4t2(l - tl)R,3R33{RnR3i + R12R32)} (B.6) 



+t/3(l + ^^)-R?3 [ipiR32 - R31) + -^31-^32(1 - 

+t^(l + t^p)Rl^ [ipiRu - Rii) + RuRui^ - 



This expression must be zero subject to the constraint flB.2p in order to have CP conser- 
vation. Thus, when ImAs 7^ and tan (3 7^ 1 we have CP conservation in the following 
cases 

• Lower mass degeneracy Mi = M2 < M3 when R^i = tan/3i?32 7^ and /x^ = 
Rl,Ml + {Rj, + Rl,)Mi. 

• Higher mass degeneracy Mi < M2 = M3 when i?ii = tan/3-Ri2 7^ and /x^ = 
Rl,Ml + {Rl, + Rl,)Ml 

• The case of nondegenerate masses Mi < M2 < M3 yields several possibilities for CP 
conservation. To see this we rewrite flB.2p as 

Mj - M| ^ R.sjRn-tanPRu) , 
Ml -Ml i?33(^3i-tan/5i?32) ^ ' ' 

which must be a positive quantity in the case of nondegenerate masses. Solving this 
equation for M| — M| and substituting into (1B.6P we end up with 

v^RUR,,-t,R,,) 5(1371) - + 

X (i?21 - tpR22) {R2I + tpR22) = 0. (B.8) 

The solutions of this equation that yield CP-conservation are those that also imply 
positive values of (M| - Ml)/{Ml - Mf)). They are 

1. i?ii + tpRi2 = 

2. i?21 + ^/3-R22 = 

3. -R31 + t/3-R32 = 

provided 

Ri3{Rii - tan/?i?i2) , , 

i?33(i?3i-tan/?i?32) ■ ^ ' ' 
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B.5 CP conservation and mass degeneracy 

Summarizing the results presented in this Appendix in terms of mass degeneracy we see 
the following: 

• When all three masses are degenerate, CP is conserved because ImA5 = 0. For 
this to happen the two input masses must be equal, Mi = M2. In most cases 
(when it!33(i?3i — i?32tan/3) 7^ 0) this will lead to M3 being equal to Mi and M2. If 
-Rssl-Rai ~ i?32tan/3) = one cannot determine M3 from the input parameters, but 
it can be arbitrarily chosen equal to the two other masses. 

• When there is only lower mass degeneracy. Mi — M2 and i?33(i?3i — it!32tan/3) = 
one cannot determine M3 from the input parameters, but it can then be arbitrarily 
chosen to be higher than the two other masses; Mi — M2 < M3. CP is then 
conserved if either: 

- i?33 = or Rsi = R32 = 0, 

- R31 — R32 7^ and tan P — 1, 

- R31 = R32 tan /3 7^ and tan /3 7^ 1 and //^ = -R33MI + {R^^ + -R32) M|. 

• When there is only higher mass degeneracy, CP is conserved in some special cases. 
For this to happen we must first choose Mi < M2 and i?i3(i?ii — i?i2tan/3) = 0. 

Then, if R^siRsi - 7?32tan/5) =^ 0, M2 = M3. If RssiR-si - i?32tan/3) = 0, one 
cannot determine M3 from the input parameters, but it can be arbitrarily chosen 
to equal M2; Mi < M2 = M3. CP is then conserved if either: 

- Ris^O or Rii = R12 = 0, 

- Rii — R12 7^ and tan (3 — 1, 

- Rii = i?i2tan/3 7^ and tan /3 7^ 1 and = Rl^M^ + (i?^^ + Rl2)Ml 

• There are also cases of CP conservation in the mass non-degenerate case Mi < 
M2 < M3. For this to happen we must first choose Mi < M2 and cki, 0:2, 0:3 and 
tan /3 in such a way that M3 > M2 or in a way such that M3 cannot be determined 
from the input parameters. Then M3 can be arbitrarily chosen to be higher than 
M2. In both cases Mi < M2 < M3. CP is then conserved if either: 

- Rii — R12 — R33 — 0, 

- R31 — R32 — R13 — 0, 

- Ri3 — R33 — 0, 

- Rii — R12 and it!3i = i?32 (excluding the three cases already mentioned above) 
and tan P — 1, 

- Rii — R12 7^ and i?33 = and tan P — 1, 

- i?3i = i?32 7^ and Ri3 — and tan P — 1, 

- Rii + R12 tan P — or R21 + R22 tan P — or R31 + R32 tan P = (excluding 
the three first cases mentioned above) and tan P ^ 1 and ji"^ takes on special 
values. 



34 



C Different basis 



For the purpose of determining the electroweak parameters T and 5* in Sec. El we need 
to relate the rotation matrix R defined by (2.3) in [20] and the O defined by (59) in 
[25]. First, we find the U and V (in order to emphasize that the inert doublet is not yet 
included, we here adopt a subscript "2HDM") of [25] in the basis adopted in [20], where 
the doublets are denoted by $i (see Eq. 



For (pi defined in the basis in which only (pi has a non-zero v.e.v. (the "Higgs basis"' 
we have 



{v + H + tG^)/V2 



iR + tI)/V2 J 



{S + tA)/V2 



■ (c.i) 



The following transformation relates and (pj: 



Cf3 Sp 



$1 
$2 



(C.2) 



with 03 = Tj. Since in ( IC.lj) the charged Higgs bosons are mass eigenstates (according to 
(15) in [25]) we obtain 

TJ _ / 
tV2HDM — 

The matrix V2HDM is defined through (see Eq. (60) in [25]) 



(C.3) 



V2 + iX2 



2HDM 



Hi 



(C.4) 



where Hi are mass eigenstates. The rotation matrix R is defined by Hi = RijTjj (see 
Eq. (12.71) ). Invoking the relation 



Xi 

X2 



C/3 -S() 
Sp 0/3 



(C.5) 



and replacing in the L.H.S. of ( ]C.4[) 771^2 and xi,2 by the mass eigenstates G^ and Hj, then 
V2HDM is identified as: 



2HDM 



icfs Rii - iSfjRis R21 - ispR2z R31 - is/sRss 
isp R12 + icpRij, R22 + icpR2'i R32 + icpRss 



(C.6) 
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